Electronic and magnetic properties of K2CUP2O7 - a model S=l/2 Heisenberg chain 

system 
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The electronic and magnetic properties of K2CUP2O7 were investigated by means of susceptibility, 
specific heat and "^^P nuclear magnetic resonance (NMR) measurements and by LDA band structure 
calculations. The temperature dependence of the NMR shift K{T) is well described by the S = ^ 
Heisenberg antiferromagnetic chain (IfAF) model with nearest neighbor exchange Ji ~ (141 ± 5) 
K. The corresponding mapping of an LDA-derived tight-binding model leads to Ji^^ ~ 196 K. 
The spin lattice relaxation rate l/Ti decreases with temperature below 300 K but becomes nearly 
temperature independent between 30 K and 2 K as theoretically expected for an 5* = 5 HAF chain. 
None of the investigated properties give any evidence for long range magnetic order above 2K, in 
agreement with the results of the band structure calculation, which yield extremely weak exchange 
to the next nearest neighbor (NNN) and a very small and frustrated inter-chain exchange. Thus, 
K2CUP2O7 seems to be a better realization of a nearest neighbor S — ^ HAF chain than the 
compounds reported so far. 

PACS numbers: 75.10.Pq, 75.40.Cx, 76.60.-k, 76.60.Cq 
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I. INTRODUCTION 

One-dimensional (ID) spin systems have attracted 
considerable attention because of their intriguing ground 
states where quantum fluctuations play a crucial role. 
Much excitement in their magnetism has been caused 
by the theoretical prediction that the integer-spin chains 
have an energy gap in the excitation spectrum,^ while 
the half-integer spin chains have a gapless excitation 
spectrum. ^'■^ In particular, ID S — ^ Heisenberg antifer- 
romagnetic (HAF) systems are interesting since enhanced 
quantum fluctuations due to reduced dimensionality and 
low spin value impede long range magnetic order (LRO). 
In the last decade, theoretical studies on these systems 
have achieved a remarkable progress but real material 
realizations for such compounds are limited to date. 

The S — ^ chains formed via direct linkage of CUO4 
units can be grouped into two categories: in one group, 
the chains are formed by edge-sharing CUO4 units while 
in the other compound family they are built from corner- 
sharing CUO4 units. CuGeOs and Li2Cu02 belong to 
the former category^, where the nearest neighbor (NN) 
interaction Ji and the next-nearest neighbor (NNN) in- 
teraction J2 have comparable strength (|j^ « 1) thereby 
leading to a strong frustration within the chain. In such 
a scenario, various ground states are possible. On the 
contrary, in Sr2Cu03 which has corner shared CUO4 
units ~ which greatly reduces the in-chain 

frustratiori^ii^. With Ji > Tat, SraCuOs is a bet- 
ter quasi-lD S = chain system than CuGeOa and 
Li2Cu02. Recently another system, Sr2Cu(P04)2 con- 
taining isolated CUO4 units (neither edge- nor corner- 
shared) was reported to have ^ « and a much 
reduced T/v of 85 mK— The presence of inter- 
chain coupling (ic) normally provides LRO, while frus- 
trating ic reduces the tendency to LRO. Along with Ji 



3> Tn, the presence of frustration between the chains 
strongly influences the ID nature of this system, making 
it a nearly perfect realization for the ID NN only HAF 
chain model. Unfortunately the sample quality was not 
very good, thereby reducing the chains to a finite length. 
The analytical solution by Bethe ansata^. provides a clear 
picture of the magnetic (thermodynamic) properties of 
a very good ID material. Thus, deviations from these 
predictions can be assigned to other degrees of freedom 
(anisotropy, Dzyaloshinskii-Moriya (DM) and other spin- 
spin interactions) allowing to analyse their influence on 
the ground state behaviour. Therefore, strong interest 
still prevails in the experimental community to find pos- 
sible ways to synthesize a system which keeps the in-chain 
geometry of Sr2Cu(P04)2, but with a smaller inter-chain 
coupling which is still frustrating, along with an improved 
sample quality. 

Formaly, K2CUP2O7 belongs to the family of 
M2CUP2O7 (M=Li, Na, K) compounds, though the ar- 
rangement of the magnetic CUO4 building blocks is mod- 
ified within the family. The magnetic properties of 
Na2CuP207 and Li2CuP207 were reported earlier.— 
Na2CuP207 and Li2CuP207 both have an in-chain ex- 
change coupling Ji of about 28 K and they undergo mag- 
netic ordering at around 5 K. Unlike the other members 
of this family, our experiments reveal that K2CUP2O7 
shows a much stronger in-chain exchange interaction 
(Ji ~ 141 K) but no ordering down to 2 K. Due to such 
a wide T/v-Ji temperature range it can be considered as 
a very good candidate on which to test the theoretical 
predictions of both static and dynamic properties of an 
ID HAF model. 

In this paper, we present susceptibility, specific heat, 
■^^P NMR results and first principles band structure cal- 
culations on K2CUP2O7 in order to shed light on the ID 
character of the system. The manuscript is organized as 
follows: In the next two sections, we describe the struc- 




FIG. 1: K2CUP2O7 crystal structure. The CUO4 square 
planes (red) share their edges with the PO4 tetrahedra (green) 
to form [Cu(P04)2]^ linear chains propagating along c- 
direction. The potassium cations reside in between the chains 
(not shown here). 



tural aspects, measuremGnt procedures and the theoreti- 
cal methodology. In section IV we present the experimen- 
tal results followed by the electronic structure analysis. 
In the discussion section, we address the issues concern- 
ing the almost ideal ID nature of this compound inferred 
from the ^^P NMR and band structure calculations. At 
the end we give a short summary and conclusions. 



II. STRUCTURE 

K2CUP2O7 crystallizes in an orthorhombic unit cell 
with space group Pbnm. The reported lattice constants 
are a = 9.509 A, b = 14.389 A, and c = 5.276 In the 
crystal structure, isolated quasi planar CUO4 units are 
linked by PO4 tetrahedra, forming [Cu(P04)2]oo chains 
propagating along the crystallographic c-direction. A 
scheme of the basic building blocks of the crystal is shown 
in Fig. [T] The super-exchange between Cu^"*" ions in 
K2CUP2O7 will be similar to that of edge sharing CUO4 
chains with every second CUO4 unit cut-off. In such a 
situation, J2 (NNN) in the edge shared system will be- 
come Ji (NN) in our system (of the order of w 100 K), 
and J2 in our system will be comparable to J4 (fourth 
neighbor interaction along the chain) in the edge shared 
system, which is known to be quite negligible. The chains 
are well separated from each other since the potassium 
cations, reside in between the chains. The inter- 

chain interactions are expected to be very weak, unlike 
Sr2Cu(P04)2 since the chains do not lie in the same 
plane. Magnetic properties of this compound have not 
been reported yet. 



III. METHODS 
A. Experimental 

Polycrystalline K2CUP2O7 was prepared by solid state 
reaction techniques using K2CO3 (99.9% pure), CuO 
(99.99% pure) and NH4H2PO4 (99.9% pure) as start- 
ing materials. The stoichiometric mixtures were fired at 
640°C for 60 hours in air, with one intermediate grinding. 
The samples were characterized using a STOE powder 
diffractometcr with a Cu target {Xav — 1.54182A). The 
powder pattern evidenced single phase material. The 
lattice parameters obtained using a least-square fit pro- 
cedure are a = 9.541(3) A, 6 = 14.407(5) A, c = 5.253(1) 
A, close to the previously reported valuesiJ^^ 

Magnetization (M) data were measured as a function 
of temperature T between 2 K and 400 K in fields up to 5 
T on powder samples in a commercial (Quantum design) 
SQUID (superconducting quantum interference device) 
magnetometer. Specific heat Cp(T) measurements were 
performed on a pressed pellet using relaxation method in 
a commercial PPMS equipment (Quantum design). NMR 
measurements were carried out using pulsed NMR tech- 
niques on ^^P nuclei (nuclear spin I = ^ and gyromag- 
netic ratio 7/27r = 17.237 MHz/Tesla) at 70 MHz which 
corresponds to an applied field of about 40.6 kOe. Spec- 
tra were obtained by Fourier transform (FT) of the NMR 
echo signal using a 7r/2 pulse with width of about 2 /is. 
The NMR shift K{T) = [v (T) - Vref] /i^ref was deter- 
mined by measuring the resonance frequency of the sam- 
ple {v (T)) with respect to a standard H3PO4 solution 
(resonance frequency Vref)- The spin-lattice relaxation 
rate (l/Ti) was measured by the saturation recovery 
method. Our attempt to see the ^^Cu signal was not 
successful due to the fast relaxation at the magnetic site. 



B. Theory 

The bandstructure calculations presented here utilized 
version 5.00-18 of the full-potential local orbital band 
structurei2ii^ (FPLO) method. The structure, lattice 
constants and atomic positions have been taken from 
experiments^. The core states have been treated fully- 
relativistically while the semi-core (K: 3s3p, Cu: 3s3p, 
P: 2s2p) and the valence states (K: 3d4s4p, Cu: 3c?4s4p, 
P: 3s3p3(i, O: 2s2p3d) are treated scalar-relativistically. 
The extent of the valence basis functions was optimized 
with respect to the total energyi^. The Brillouin zone 
sampling was based on 216 fc-points in the irreducible 
wedge. We have used the Perdew and Wang^S flavor of 
the exchange and correlation potential when performing 
the calculations within the local density approximation 
(LDA). 
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IV. RESULTS 
A. Susceptibility and specific heat 

Bulk magnetic susceptibility x{T) {— M/H) was mea- 
sured as a function of temperature in an applied field 
of 5 kOe (Fig. [J). As shown in the figure, the sam- 
ple exhibits a shoulder at about 100 K, suggesting a 
maximum in this temperature range, which is a hall- 
mark of low-dimensional magnetic interactions. With 
further decreasing temperature, xi^) increases again in 
a Curie- Weiss manner. Common sources for such an in- 
crease are the presence of paramagnetic impurities or fi- 
nite length chains due to disorder. However, an intrinsic 
Curie like tail is expected in a quasi ID 5 = ^ chain, 
if DM interaction and/or a staggered factor anisotropy 
are present.^^''^^ Among the 3d systems, Cu benzoate^^ 
and pyrimidine Cu dinitrate^ are well known examples. 
No obvious features associated with LRO were observed 
for 1.8 K < T < 400 K. 

In order to fit the bulk susceptibility data, we decom- 
posed X into 



X = X0 + 7f7-^ ^Xspin{T) (1) 

' "imp 

where, the first term xo is temperature independent and 
accounts for the diamagnetism of the core electron shells 
and Van-Vleck paramagnetism of the open shells of the 
Cu^+ ions. The second term S.'^T'^ is the low T Curie- 
Weiss contribution due to paramagnetic species in the 
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FIG. 2: Magnetic susceptibility [M/H) vs. temperature T 
for K2CUP2O7 (open circles) in an applied field of 5 kOe. The 
solid line is best fit of the data to Eq. 1. Spin susceptibility 
Xspin is plotted (open triangles) after subtracting the Curie 
contribution. In the inset, Cp/T is plotted as a function of 
for the low temperatures. Solid line represents a linear fit. 



sample. Xspin{T) is the spin susceptibility for a uniform 
5 = i ID HAF system, which is known quite precisely 
over the whole measured temperature range. We took 
the Johnston's expressions^ valid for 5 x 10^^^ < 7^ < 5. 



A fit of Eq. 1 to the experimental data leads to the 
parameters xo — 2.7x10^'* emu/mole, Cimp — 0.026 
emu.K/mole, 9imp — 0.7 K, and Ji ~ 145 K. Because 
of large Curie tail, such a fit is unstable at high tem- 
peratures. Therefore we reduced the number of fitting 
parameters by fixing g = 2.2 (obtained from NMR shift 
analysis)22i. The value of xo is comparable to that found 
in Sr2Cu03i. The Curie contribution present in the sam- 
ple would correspond to a defect spin concentration of 6.9 
% assuming defect spin S = Alternatively, it would 
correspond to the contribution expected for finite chains 
with a length of « 15 spins according to the calcula- 
tion of Schmidt et al.r^. We add in Fig. 2 a plot of 
Xspin as a function of temperature, after subtracting the 
Curie contribution. Now a pronounced broad maximum 
around 100 K is well resolved. The small peak at 50 K 
corresponds to a parasitic contribution of adsorbed oxy- 
gen . Unfortunately, we failed to improve the quality of 
the sample, since further annealing (in flowing Ar or vac- 
uum) did not lead to any noticeable change in the XRD 
pattern. 



The information about the lattice dimensionality can 
be obtained from the low temperature specific heat 
Cp{T). For a S — ^ chain, one expects a linear 
term, whereas for a square lattice, one expects a lead- 
ing quadratic term in the absence of excitation gap. In 
the Cp/T vs. plot (shown in the inset of Fig. [2]) in 
the temperature range 5K<T<12K, one can see 
clearly that the plot follows a straight line over a con- 
siderable temperature range. This indicates that Cp(T) 
is a sum of a linear and a cubic contribution. Since the 
cubic term corresponds to the expected contribution of 
the phonons, this demonstrate that the leading term of 
the magnetic contribution is linear in T . For an 5 = ^ 
HAF chain, theoretical calculations^!^ predict for low 
temperatures (T < 0.2Ji); ^ ^ ^ = jtheo- With 
Ji ~ 141 K obtained from the NMR shift K{T) analysis 
(presented later), this corresponds to a value ■^theo — 0.04 
J/K^ mole. Fitting the measured data in the range 5 K 
< r < 12 K, we obtained 7oxp ^ 0.042 J/K^mole. The 
value of 7exp is close to the value predicted theoretically. 
This strongly supports the quasi ID nature of the spin 
system in K2CUP2O7. No other anomaly was observed 
down to 2 K in Cp[T), suggesting the absence of mag- 
netic order. Below 5 K, our experimental data deviates 
upward from the straight line fit (inset Fig. [2]). This is 
likely related to the contributions of paramagnetic impu- 
rities and of chain ends, which should follow a C « A/T"^ 
behavior. 
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FIG. 3: -^^P NMR spectra at different temperatures T. Solid 
line represents the nonmagnetic "^^P reference. 



B. ^^P NMR 

Although our analysis of the susceptibility and specific 
heat suggests the presence of a 5' = ^ ID HAF system, 
the evidence is only of preliminary nature because of the 
large Curie tail in x(2^) &t low T. In order to gain a more 
reliable insight into this system, we turned our attention 
towards NMR results. NMR has the advantage to be 
much less sensitive to contribution of defects or impuri- 
ties, because usually only the nucleus on an undistorted 
site contributes to the narrow NMR line. As shown in the 
crystal structures (Fig. [J), K2CUP2O7 has two incquiv- 
alent ^^P sites which are coupled inductively to Cu^+ 
ions in the chain. Therefore ^^P NMR can probe accu- 
rately the low-lying excitations in the spin chain. Our 
•^^P NMR spectra consist of a single spectral line as is 
expected for / = ^ nuclei (Fig. [3]). Although there are 
two inequivalent ^^P sites present the crystal structure, 
a single resonance line implies that both the '^^P sites 
in this compound are nearly identical. With decreasing 
temperature, the NMR line shifts away from the Larmor 
frequency, broadens paramagnetically but the overall line 
shape remain same down to 2 K. The asymmetric shape 
of the spectra corresponds to a powder pattern due to an 
asymmetric hyperfine coupling. 

The temperature dependence of the NMR shift K is 
shown in Fig. [H With decreasing temperature K{T) 
increases paramagnetically, then passes through a broad 
maximum at 110 K, which is an indicative of short-range 
correlations, and decreases again smoothly towards low 
temperatures. Below -j^ ~ 0.028, K{T) shows a much 
steeper decrease towards zero. As mentioned before the 
NMR has an advantage over bulk susceptibility. One 
measures accurately the Xspin by NMR shift without suf- 
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FIG. 4: Temperature dependence of the ^^P NMR shift K{T) 
of K2CUP2O7. The solid line is the fit with Eq. 2 in the tem- 
perature range, 9 K < T < 300 K and further extrapolated 
down to K. Lower inset shows K vs. T on a logarithmic tem- 
perature scale for improved visualisation of the low-T data. 
In the upper inset, K vs. x is plotted with temperature as an 
implicit parameter and the solid line is the linear fit. 

fering from the contribution from the free spins and ex- 
trinsic foreign phases, which limits the accuracy of the 
bulk susceptibility measurements. Therefore, it is more 
reliable to extract the magnetic parameters from the tem- 
perature dependence of the NMR shift rather than from 
the bulk susceptibility. The conventional scheme of cor- 
relating K{T) and x{T) is to plot, K vs. Xspm with T 
as an implicit parameter. Then the slope yields the av- 
erage hyperfine coupling A^f between the ^^P nucleus 
and the two nearest-neighbor Cu^+ ions. In the case of 
K2CUP2O7, because of the low T Curie-tail in x(T) this 
K{T) vs x{T) plot shows a straight hue only for T > 
110 K (see upper inset of Fig. U) Nevertheless, we can 
estimate A^f ~ (4400±400) Oc/fis, which is about two 
times stronger than in other phosphate systems^^'iS.. We 
then determine Ji and g by fitting the temperature de- 
pendence of K to the following equation, 

K = Ko+(^-^^Xsp^n{T,Jl) (2) 

where Kq is the temperature independent chemical shift. 
As shown in Fig. [H the K{T) data fit rather well to Eq. 
[5] in the temperature range 9 K < T < 300 K. Using 
Afif — 4400 Oe/ jiB (obtained from the vs x analysis) 
we obtained Kq ~ — SOOpprn, Ji ~ (141 ± h)K , and g ~ 
2.2. Below about 8 K, K{T) show a significant deviation 
from the fit (see lower inset of Fig. 2]) which shall be 
discussed later. 

The temperature dependence of '^^P l/Ti is presented 
in Fig. [51 For a / = ^ nucleus the recovery of the 
longitudinal magnetization is expected to follow a single 



4 



exponential behavior. In the experiment, we indeed ob- 
served single exponential behaviour down to 30 K while 
below 30 K it fitted nicely to the stretch exponential with 
a reduced exponent. We didn't observe any anomaly or 
divergence in l/Ti{T) down to 2 K which indicates the 
absence of magnetic ordering. For 2 K < T < 30 K, 1/Ti 
almost remains constant while for T > 30 K, it increases 
strongly with temperature. A slight shoulder is visible in 
the plot 1/Ti versus InT around 180 K, i.e just above Ji, 
suggesting a regime crossover in this temperature range. 



C. First Principles and Tight Binding 




Collected in Fig. [6] are the non-magnetic band struc- 
ture and the density of states (DOS). The antibond- 
ing band is made from a half-filled Cu - 3^3.2 _j,2 - O - 
2pcr molecular orbital belonging to the CUO4 plaquettes. 
The band structure shows strong dispersion of about 55 
meV parallel to the chain direction TZ, XA, but is nearly 
dispersion-less within the crystallograhic a — b plane {i.e. 
perpendicular to the direction of the chains), indicative 
of the strong ID character in this system. Within LDA, 
we get a metallic behavior, though the system is an insu- 
lator, suggested by its blue colorili This is a well known 
artifact of LDA wherein the effect of strong correlations 
is under estimated. Including the strong correlations for 
the Cu 3c? states will open up the insulating gap. This 
can be achieved via two possible ways: (a) performing 
an LDA-|-f7 calculation self-consistently, (b) mapping the 
results from LDA first to a tight-binding model (TBM). 
The TBM is then mapped on to a Hubbard model, and 
subsequently to a Heisenberg model because the system 
belongs to the strong correlation limit C/ 3> t (t is the 
leading transfer integral) at half filling. 



0.5 

0.4 

^ 0.3 

'w 

E 

0.2 
0.1 
0.0 



40 D 

3 



in o o P oopc" 



0° 



50 100 150 200 250 

T(K) O 



CD°0 § 000 C0° 



10 



T(K) 



100 



FIG. 5: Spin- lattice relaxation rate 1/Ti vs. temperature T 
measured at 70 MHz. In the inset, 1/{KT\T) is shown all 
over the temperature range. 



FIG. 6: (color online) Left panel: The non-magnetic band 
structure of K2CUP2O7 along the high-symmetry directions of 
a standard orthorhombic unit cell. X—A and Z—V denote the 
directions along the chain, (i.e. along crystallographic c axis) 
wherein the dispersion is the largest while the perpendicular- 
to-chain directions V — X and T — Y are rather non-dispersive. 
There are four copper sites per unit cell leading to four half- 
filled metallic bands at the Fermi level, well separated from 
all other bands. Right panel: The total DOS for K2CUP2O7 
along with the Cu — d^2_y2 orbitally resolved contributions. 
Comparison with the band structure clearly elucidates that 
the half-filled metallic band at the Fermi level is primarily 
from Cu — d^2_y2 and O — pa molecular plaquette orbital. 



In order to better understand the microscopic mag- 
netic interactions in this system, we have followed 
the second option mentioned above. Firstly we have 
considered the TBM, 



H = 



E 



(3) 



where ty are the hopping integrals, Cj^Cj,^ are anni- 
hilation and creation operators. From the LDA band 
structure, we have extracted only the 4 antibonding Cu 



3d;. 



bands and performed a fit to the TBM. All 



the hopping paths considered in our model are shown 
in Fig. [71 The hopping integrals were calculated using 
the steepest descent method. The resulting parameters 
which provided the best fit to the LDA band structure 
(Fig. [g are collected in Table U The strength of the NN 
hopping ti along the chain is two orders of magnitude 
larger than all the other hoppings, attributing the strong 
one-dimensionality to the interaction of the plaquettes 
along the chain. The individual exchange constants are 
calculated using J^^^ = AtjjUeff. The value of Ueff 
for K2CUP2O7 is set to 4.5 eV, same as the choice in 
related ID compounds (Sr,Ba)2Cu(P04)2^. The total 
exchange constant is given by J*°*°' = jAfm _|_ jFM _ 
The J values collected in Table U are indicative of the 
jAFM Q^iy^ jj^ K2CUP2O7 the CUO4 plaquettes are sep- 
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arated from each other, with no corner- or edge-sharing 
oxygens. The absence of direct connections between the 
plaquettes largely suppresses ferromagnetic interactions 
between the copper sites. The strength of the FM inter- 
actions was shown to be very small in the related systems 
(Sr,Ba)2Cu(P04)2^^, and this result should equally hold 
for our system. 




FIG. 7; (color online) The various hopping paths considered 
in our tight-binding model to reproduce the half-filled metal- 
lic LDA band are shown here. We have considered two hop- 
pings along the chain (ti,t2) and four inter-chain hoppings 
{tici,tic2,tic3,tic4) in all. Starting from one chain, Cu in the 
next neighbor chain are shifted by half of the in-chain Cu-Cu 
distance, leading to two identical interactions to two Cu-spins 
in each adjacent chains. 
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TABLE L Hopping parameters (in meV) and the correspond- 
ing exchange constants J ( in K) from an effective one-band 
tight-binding model. The hopping paths used are indicated 
in Fig. [3 
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FIG. 8: (color online) The superposition of the total band 
structure from the FPLO density functional theory (DFT) 
calculations along with the calculated eigenvalues from the 
four-site one-band TBM. The TBM fit is quite consistent with 
the DFT resuhs. 



V. DISCUSSION 

The quality of the fit for the NMR shift supports the 
presence of an 5 = ^ HAF chain. The exchange cou- 
pling ( Ji ~ 141 K) is comparable to that found in other 
phosphates^ and in nice agreement with the effective NN 
super-exchange constant, — 196 K obtained from 

the TBM. The slight over-estimation may stem from the 
fact that Uef / is not exactly known, along with some FM 
contributions. Such a over-estimation of J by LDA is 
well known and is also observed in (Sr,Ba)2Cu(P04)2^ 
The next nearest neighbor (NNN) super-exchange con- 
stant ^2^"^ — 0.04 K, is extremely small, so frustration 
coupling, if any should be negligible. The ratio of first 
and second neighbor in-chain coupling is jf ^"^/jf"^"* > 
5000, the largest found in 5 = i chain system. The ra- 
tio of in-chain to the strongest (frustrating) inter-chain 
couphng is Jf J,^f ^ > 1000, which is two orders of 
magnitude larger than in Sr2Cu(P04)2 {Jt^^/Jt^^ ^ 
70)^^ and one order of magnitude larger than Sr2Cu03 
{Jt^^/Jt^^ 500)^ and making K2CUP2O7 an 
even better realization of ID HAF behavior than the 
Sr analogue. We have estimated the Neel temperature 
of K2CUP2O7 (T^) by adapting a simple mean field 
approximation^ and comparing it to the Sr analogue 
(T^' ). Assuming that the anisotropy are the same in 
both compounds we can write according to Ref. 

rpK I jK jK 

where, Tj^'" « 85 mKiS. Since the incompleteness in the 
mapping should be the same for both the compounds, 
we can directly compare the values of T^r. We have 
used Jf'' = 187 K, jfj' = 0.23 K (geometrical average 
of the 2 inter-chain exchanges), = 196 K and = 
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J' 



0.037K. Plugging these val- 



ues in the above equation we get 37 mK, which is 

a factor of two smaller than the Sr analogue. This value 
sets a upper bound, because the value of Jic is calcu- 
lated using a mean-field approximation. The absence of 
anomalies in K{T), l/Ti{T) and the invariant spectral 
shape down to low temperatures rules out the possibil- 
ity of LRO down to 2 K, consistent with the theoretical 
prediction of an extremely low T/v. In K2CUP2O7, each 
Cu04 plaquettes in one chain has two identical neighbors 
in each adjacent chains, with the same exchange interac- 
tions. For a strong in-chain AFM exchange, this leads 
to a complete frustration of the inter-chain interactions. 
Quantum fluctuations are therefore enhanced here which 
suppresses Tpf to lower values. 



For a ID S' 



HAF, theoretical calculations pre- 



dict a weak logarithmic decrease of xiT) upon approach- 
ing T = K (see Fig. S]). The lower inset of Fig. g] 
shows clearly that the decrease we observe in K{T) of 
K2CUP2O7 is much more pronounced . The susceptibil- 
ity of a ID 5 = i HAF at T = is exactly knowi 



X{T = 0) = 



Then, K{T) at zero temperature 
0) - i^o 



(^1) 

determined 



can be written as Ktheo{T 

Using the parameters {Kq, Aj^f, g, and Ji 
from our K{T) analysis, Ktheo{T — 0) was calculated 
to be 260 ppm. However, our experimental value at 2 
K is much lower, only 180 ppm, and K{T) is still de- 
creasing steeply with T at this temperature. Thus, at 
the quantitative level, the theoretically predicted loga- 
rithmic term fails to describe our experimentaly observed 
decrease. A similar feature has been found experimen- 
tally in NMR on SrzCuOa below T/Ji ~ 0.015^0 and 
3ip NMR on (Sr,Ba)2Cu(P04)2 below T/Ji ~ 0.003.^ In 
both the cases, it is argued that the decrease is unrelated 
to the onset magnetic order or a spin-Peierls transition. 
In K2CUP2O7, we neither saw any indication of order- 
ing in l/Ti(T) nor in Cp(T). Similarly there is no sig- 
nature of exponential decrease (or singlet ground state) 
observed in 1/Ti (T) as is expected for a spin-Peierls tran- 
sition. One possibility to account for this drop in K[T) is 
the presence of DM interaction arising from the fact that 
in K2CUP2O7 there is no center of inversion symmetry- 
relating two neighboring copper atoms along the chain. 
In the presence of DM interaction, application of a mag- 
netic field parallel to the DM vector opens a gap A in 
the magnetic excitation spectra. Since in a first approxi- 
mation A increases with B to the power 2/3, it shall not 
be visible in the B = Q specific heat data or the low field 
susceptibility shown in Fig. [2l However, further experi- 
ments are needed to confirm or discard this explanation. 

In NMR valuable information on the dynamic of low- 
energetic spin excitations can be gained from the analy- 
sis of the temperature dependence of nuclear spin-lattice 
relaxation rates. Therefore, it is essential to analyze 
l/Ti{T) carefully which yields information about the 
imaginary part of the dynamic susceptibility xi^-i^)- 
Thus 1/Ti should include contributions from both the 



uniform [q = 0) and staggered {q = spin fluctua- 
tions. A theoretical analysis by Sachdev^ shows that 
the staggered component is dominant at low tempera- 
tures [T <^ Ji). Indeed, the uniform component leads 
to 1/Ti oc T, while the staggered component gives l/Ti 
= constant. Monte Carlo calculations by Sandvik sup- 
ported the validity of these results over an appropriate 
temperature range. '^■^ As shown in Fig. [5l our exper- 
imentally observed constant l/Ti at low temperatures 
{2K < T < 30K) suggests the dominance of staggered 
fluctuations at low temperatures. The '^"'^P form factor 
for such systems is defined in Ref.[^. Since ^^P is located 
symmetrically between the Cu ions, the fluctuations are 
expected to be filtered out provided the hyperfine cou- 
plings are equal. Moreover in this case we still have a 
significant remnant contribution from q = ±7r/a which 
plays a dominant role at low temperatures. The possible 
origin of the remnant staggered fluctuations could be the 
unequal hyperfine couplings. In fact, such features have 
been previously observed in a few other ID S = ^ HAF 
systems where Ji ^ T/v i^'^*^'^'^ The strong increase 
of 1/Ti with T above 30 K could be attributed to the 
uniform fluctuations. When the dominant contribution 
is from q — 0, then one expects a constant l/{KTiT). 
In K2CUP2O7, for T > 30 K, we rather observe a weak 
temperature dependency (inset of Fig. [5]) which might be 
due to some remanent contributions from the staggered 
fluctuations and/or some additional relaxation mecha- 
nisms. This weak temperature dependence as well as the 
reduction of the slope in the plot 1/Ti versus T above 
180 K might be related to spin diffusion as observed in 
a-VO(P03)2 above Ji^?. 

The relaxation rate due to staggered fluctuations can 
be calculated following the prescription of Barzykin— . 
For the purpose of comparison of theory with experi- 
ment he defined the normalized dimensionless NMR spin- 
lattice relaxation rate at low-temperature (l/rL)^^^,,^ = 
.^•^^ « 0.3, where Ath is Ahj (2/i7/27r). Assuming the 
fiuctuations to be correlated because of the exchange Ji 
along the chains, 1/Ti can be written as l/Ti = -t-tta^- 

Using this expression, (1/Ti) at the '^^P site was calcu- 
lated to be about 129 sec,^^ whereas our experimental 
value is 50 sec,-^ in the 2 K < T < 30 K range. The 
experimental value is about two times smaller than the 
theoretical value, likely due to the effect of the geometri- 
cal form factor. Further on a weak logarithmic increase 
in 1/Ti is theoretically expected at low temperaturesi^ 
In the present case, our measurements were not done 
down to low enough temperatures to rise this contribu- 
tion above our experimental error. 



VI. CONCLUSION 

In conclusion, our experimental and theoretical studies 
of K2CUP2O7 demonstrate that this compound presents 
a very uniform and strongly ID S — HAF chain system. 
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Our NMR, specific heat and susceptibility measurements 
show good agreement with theoretical predictions for a 
ID 5 = ^ HAF chain. Thus, the temperature depen- 
dence of the NMR shift K can be fitted in good agree- 
ment with model calculations giving Ji ~ (141±5) K. A 
calculation of the exchange interactions using a TBM fit- 
ted to the results of ah initio LDA band structure calcu- 
lations, subsequently mapped onto a Heisenberg model, 
leads to a quite similar value of Ji ~ 190 K and evidences 
extremely weak NNN interactions as well as extremely 
weak and "frustrated" inter-chain exchange of the order 
of J"" ^ 0.1 K, in contrast to Sr2Cu(P04)2 with J"-" ~ 9 
Kii and to SrsCuOa "non-frustrated" J'" ~ 10 Kl. The 
absence of any evidence for magnetic order above 2 K in 
the experimental data confirms the weakness of the inter- 
chain exchanges. Using the TBM results, a rough esti- 
mate of the Neel temperature results in a value of Tn ^ 
37 mK, only. K2CUP2O7 has the smallest in-chain as 
well as frustrating smallest inter-chain exchanges which 
do not suppress quantum fluctuations, thereby making 
this system an even better example of a ID 5 = ^ HAF 



system than the compounds reported so far. 

We further studied the magnetic fluctuations by 
analysing the temperature dependence of l/Ti. At low 
temperatures l/Ti remains constant, in reasonable agree- 
ment with Sachdev's predictions for ID 5 = ^ HAF sys- 
tem where relaxation is dominated hy q — ±7r/a at low 
T. Below 5 K, we obtained a decrease of the NMR shift 
K which is more pronounced than that expected for a 
ID = i HAF. It's origin is not yet clear, but might be 
due to DM interactions. 
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